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Abstract 

Based on an effective Lagrangian which combines chirai SU(3) dynam- 
ics with vector meson dominance, we have developed a model for the 
forward vector meson-nucleon scattering amplitudes. We use this as 
an input to calculate the low energy part of the current-current cor- 
relation function in nuclear matter. For the isovector channel we find 
a significant enhancement of the in-medium spectral density below the 
p resonance while the p meson mass itself changes only slightly. The 
situation is different in the isoscalar channel, where the mass and peak 
position of the u meson move downward while its width increases mod- 
erately. For the 4> meson we find almost no mass shift, but the width 
of the peak increases significantly. We use these spectra as "left hand 
side" of in medium QCD sum rules and observe a remarkable degree 
of consistency with the operator product expansion "right hand" side 
in all three channels. We point out, however, that these results can- 
not simply be interpreted, as commonly done, in terms of a universal 
rescaling of vector meson masses in matter. We also compare the re- 
sulting thermal dilepton rates of our model with CERES data. We find 
satisfactory agreement but we point out that the dilepton rates from 
completely uncorrelated quark-antiquark pairs also fit the data quite 
well. Considering the strongly attractive potential for the u> together 
with its comparably small width, we discuss the possible formation of 
bound u) states in nuclei. Using a Green's function approach we inves- 
tigate the possibility of detecting these bound states in the A(d, 3 He) 
pick-up reaction. 



1 Introduction 

It is well known that the chirai symmetry of QCD in spontaneously broken. 
This leads to the octet of light pseudoscalar mesons which are the Goldstone 
bosons of this broken symmetry, and to a mass splitting between chirai part- 
ners of all hadrons. An example is the vector meson whose chirai partner, 
the axial vector meson, is twice as heavy. There are, however, strong hints 
from the quark condensate in medium QCD lattice simulations 0, [3| or 
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schematic chiral models such as the Nambu-Jona-Lasinio model ||, that this 
symmetry is expected to be restored at high temperatures and high densities. 
In order to find some signals for this chiral restoration it was proposed to study 
dileptons emitted in heavy ion collisions since they can carry information of 
the inner hot and dense zone of such collisions. Such experiments with high 
statistics and resolution are planned by the HADES collaboration at GSI. For 
high temperature there are already data from the CERES || and HELIOS 
U collaborations at CERN which observe strongly enhanced dilepton yields 
below the rho meson resonance. This was often interpreted as an indication 
of dropping vector meson masses which in addition should be a signal for chi- 
ral restoration. The Brown-Rho scaling scenario [7| and considerations based 
on bag models || support this idea and find a strong reduction of the vector 
mesons masses. Also the QCD sum rule analysis |§, which however only used 
a delta function caricature of the true spectrum, seemed to confirm these re- 
sults. From the theoretical point of view chiral restoration, however, does not 
demand a drastic reduction of the vector mesons masses in medium. Hadronic 



calculations [10], [TTJ, [ll| suggest only marginal changes of the in-medium rho 
meson mass, but a strongly increased decay width instead. It was shown 
[0, |13[ that this is also in agreement with QCD sum rule analysis if one re- 
places the delta function by the full spectrum. Taking also p-wave interactions 
into account this leads to an agreement with the CERES dilepton data. 



In a thermal model the dilepton yields of heavy ion collisions are governed 
by the thermal dilepton rates given by 

'^Imntg, p,T)f B (q,T), (1) 
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with a = e 2 /47r = 1/137. We denote the four-momentum of the dilepton 
pair by q, the temperature and the density of the system with T and p. The 
space-time volume element which radiates the dilepton pair is d A x. Apart 
from the trivial Boltzmann factor fg, only the imaginary part of the current- 
current correlation function H(q, p, T) , which we will investigate in this talk, 
determines directly the dilepton rates. 

In section 2 we review the properties of the current- current correlation 
function in the vacuum. The corresponding correlation functions in baryonic 
matter will be developed in section 3. To leading order in the density this 
requires a detailed calculation of the vector meson-nucleon scattering ampli- 
tudes. The in-medium spectra and their comparison with QCD sum rules will 
then be presented in chapter 4. In section 5 we compare the resulting dilepton 
rates with the CERES data. We also show the result from the correlation 
function which originates from uncorrelated quark-antiquark pairs. We will 
see that both scenarios agree with present data. We therefore suggest to study 
vector mesons in additional experiments under more moderate conditions, such 
as the A(d, 3 He) pick-up reaction, to gain further informations. 



2 Vector mesons in the vacuum 



Before we start looking into the in-medium properties of neutral vector mesons 
it is useful to give a brief reminder of their vacuum properties. Since the vector 
mesons are not stable they can only be seen as resonances in the current-current 
(CC) correlation function, 

= i J d 4 xe^(0\T Jfl (x)U0)\0), (2) 

where T denotes the time-ordered product and is the electromagnetic cur- 
rent. It can be decomposed as 

u ./;; • • A (3) 

into vector currents specified by their quark content: 

ft = ^(ui^u-d^d), (4) 
ft = ^{ui^u + d-f^d), (5) 

ft = -\(^)- (6) 

Current conservation implies a transverse tensor structure 

IW?) = - qj f) n(g 2 ), (7) 
with the scalar CC correlation function 

I% 2 ) = \g^M). (8) 



The imaginary part of the correlation function is proportional to the cross 
section for e + e~ — > hadrons: 

, a(e + e~ — > hadrons) 127r T rr . . . . 

R(s) = - a t^ — r = Imn s 9 

where a(e + e~ — > ^ + ^) = 47ro; 2 /3s with a = e 2 /4ir = 1/137. The vector 
mesons can be distinguished by looking at different hadronic channels with 
corresponding flavour (isospin) quantum numbers. We show for example the 
isovector current describing the p meson. G-parity demands that it decays 
only into even numbers of pions (see data in fig. la). We observe two different 
energy regions: the low energy region of the CC correlation function is very well 




Figure 1: Spectrum of the isovector current- current correlation function in 
the vacuum and the pion form factor . The solid lines shows result from our 
vector meson dominance calculation 
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described by a hadronic language through vector meson dominance (VMD). 

V) 



The extended VMD model [15] gives 
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which is shown by the solid line in Fig la. It is determined by the imaginary 



part of the vacuum p-meson self energy n™ c [To| , which comes from the decay 



of the p into two pions, and the pion form factor 
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shown in fig lb explicitly. Here ra p is the bare p meson mass and g (g p ) is the 
coupling of the p meson to the pion (photon). The ratio of g to g p is identical to 
unity in case of "complete" VMD in which all of the hadronic electromagnetic 
interaction is transmitted through vector mesons. Note in passing that the 
denominator of the pion form factor is the full p meson propagator, which 
originates from the intermediate p-meson as depicted in fig. 2 and gives rise 
to the peak in the pion form factor at the rho meson mass. Isospin violating 
processes are small but visible as pu mixing corrections in the peak of the pion 
formfactor (see Fig lb). 

In the high energy region on the other hand the measured correlation func- 
tion approaches the asymptotic plateau predicted by perturbative QCD: 
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Figure 2: Diagrammatic resolution of the pion form factor 

The threshold so lies at about 1.2 GeV 2 . This continuum stems from a picture 
where the photon couples directly to a quark-antiquark pair, which travels 
freely over a short distance and then hadronizes into the multi pion continuum. 
Note that a system of uncorrelated quark and gluons (" Quark-Gluon-Plasma" ) 
would also lead to such a plateau in the spectrum. The threshold would then 
start already at twice the current quark mass. 



3 Current correlation functions in baryonic matter 

The CC-correlation function in medium at temperature T = is defined as 



-hoc 
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where we have replaced the vacuum by \p), the ground state of infinite, isotropic 
and isospin symmetric nuclear matter with density p. We assume matter as 
a whole to be at rest. This specifies the Lorentz frame that we will use in 
the following. For q = 0, the case where the vector meson is at rest, only the 
transverse tensor structure survives and IT 00 = IT - 7 = IT = 0. We write 

Uij(u, q = 0;p) = -SijU(uj,q^ 0; p). (14) 

As a first step we apply the low density theorem [|l(| to the self energy of 
the vector mesons: 



uV(w, q = 0; p) = Iin^) - pT VN (u) + ... , (15) 
where T VN is the vector meson-nucleon scattering amplitude taken at q = 0. 



We can extend eq.(10) to the in-medium correlation function in terms of VMD 
and write IT2] 
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Figure 3: Dominant diagrams for the Compton amplitude T(uj) and the vector 
meson-nucleon amplitudes T vn {uj). In order to get from T to T VN , replace 
the (wavy) photon line by the relevant vector mesons. Diagrams (e) and (f) 
operate, in particular, in the uj meson channel. 



We wrote down here the general form for all three channels replacing g by 
gy and denoting the ratio by ay. The vacuum vector meson self energy was 
replaced by eq.(|Tj|) which includes the complete set of medium modifications. 
In order to determine the in-medium correlation function we are left with 
the task to calculate the vector meson-nucleon scattering amplitudes. We use 
an effective Lagrangian which combines chiral SU(3) dynamics with VMD. 
This Lagrangian has been applied successfully in the vacuum |TjJ and has 
been extended to incorporate meson-baryon interactions. For the baryons we 
include the complete baryon octet and decuplet (nucleons, hyperons, A's, etc.). 

The most important processes contributing to the scattering amplitudes 
are shown in fig. 3 for the p, uj and meson respectively. For the p ((f)) meson 
we draw the diagrams which survive in the limit of large baryon mass (fig. 
3a-d). The last one (fig. 3d) only contributes to the real part of the pN 
(4>N) scattering amplitude. While for the p-meson the nN and ttA loops 
govern the scattering amplitude, KY,, KA and iTE(1385) loops dominate for 
the 4> meson. We also include s-wave scattering of the kaons (fig. 3h) p"Tf . 
For the uoN scattering amplitude we only show the dominant contributions 
to the imaginary part (figs. 3e, f, g). We evaluate the imaginary parts of 
the amplitudes by cutting the diagrams in all possible ways. The real parts 
are then determined by using a once subtracted dispersion relation, with the 
subtraction constant fixed by the Thomson limit. Evaluating those diagrams 
and using eq. ( [TBI) we plot the spectra of the correlation functions for various 
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Figure 4: Calculated spectra of current-current correlation functions. The 
dashed lines show the vacuum spectra in the p, oj and channels normal- 
ized such that they can be compared directly with the corresponding e + e~ — > 
hadrons data. The long dashed and solid lines show the spectral functions in 
nuclear matter at densities p /2 and p — 0.17 fm~ 3 , as discussed in section 3. 

densities as shown by the curves in fig. 4a, b and c. We observe important 
differences between the various channels. The p meson mass decreases only 
slightly while the width increases very strongly. This causes the peak to shift 
downwards and broaden. We also see strong threshold contributions starting 
at the pion mass. The p-meson therefore dissolves in nuclear matter; it is not 
a good quasi-particle any longer. On the other hand the mass of the u meson 
decreases significantly. Its width also increases but not as strongly as for the p 
meson and so it can still be regarded as a good-quasi particle. For the meson 
there is almost no change in the peak position while the width increases. 



4 Comparison with QCD sum rules 



Up to now we only looked at the time like region of the correlation function 
where it has an imaginary part and is therefore experimentally accessible. Now 
we want to concentrate on the spacelike region. For momentum \q 2 \ 3> 1 GeV 
the correlation function is governed by perturbative QCD which and is well 
under control. The interesting region is the space like region at around Q 2 = 
—q 2 ~ 1 GeV where non-perturbative effect start to become important. The 
basic idea of QCD sum rules is to compare two ways of determining n(g 2 ) 
in this region. One way is to use a twice subtracted dispersion relation for 
each one of the channels V = p, u, <fi: 



U(q 2 ) = 11(0) + cq 2 + — I ds 

7T J 



ImU(s) 



s 2 (s — q 2 — ie) 



(17) 



where one uses the measured or calculated spectrum Imll in the time-like 
region as input. The other one is to calculate the correlation function using 
the operator product expansion (OPE): 
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Here d agrees with the value of the perturbative continuum (e.g. d=3/2 for the 
isovector channel) and the coefficients ^2,3 incorporate the non-perturbative 
parts coming from the condensates such as the gluon condensate in C2 or 
the four-quark condensate in C3. In medium these coefficients become den- 
sity dependent because the condensates change in matter and new conden- 



0.04 + 0.018(p/p ) GeV 4 and 



sates arise. We use the values 

— 0.07+0. 036(p/ p )GeV 6 proposed by Hatsuda et al. || for the p and u meson. 
They neglected higher twist condensates which are hardly known but might be 
important. For the four-quark condensates we assume as in ref.0 that ground 
state saturation holds to the same extent as in the vacuum. The coefficient 
ci is proportional to the squared quark mass; due to the small mass of the up 
and down quark we can safely neglect those contributions and set c%' u = 0. 
The situation is different in case of the meson. Because of the large strange 
quark mass Ci is no longer negligible and we take the values cf = —0.07 GeV 2 , 
4 = -0.1 + 0.01(p/p ) GeV 4 and ct = -0.07 + 0.006(p/p )GeV 6 . 

A Borel transform is used in order to improve the convergence of the OPE 
series. Comparing eq. dl7|) and (|18D after Borel transformation we end up with 
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where the Borel mass parameter should be chosen in the range -M>1 GeV in 
order to ensure convergence of the OPE side of eq. (|T9D . The "left side" comes 
from the dispersion relation of eq.([H]) and R represents the ratio (^), but 
now specified for each individual flavour channel with V — p, u, <fi. The "right 
side" is determined by the OPE. In fig. 3a, b and c we show the comparison 
between the "left side" (solid line) and the "right side" (dashed line) for the 
various channels. The consistency between "left" and "right" sides in the 
vacuum is evidently quite satisfactory. Only for a Borel mass below 0.8 GeV 
higher order terms in the OPE become important and the convergence fails. 
The agreement at normal nuclear matter density is not as excellent as in the 
vacuum but still very impressive. We therefore conclude that QCD sum rules 
and our hadronic model of the in-medium current-current correlation function 
are mutually compatible. We also point out that assuming a simple pole ansatz 
for the spectrum at finite densities can lead to erroneous interpretation of the 
in-medium masses. This is very obvious in the isovector channel. While our 
model gives only a small change of the rho mass a simple pole ansatz leads to 
a reduction of m p by more than 10 percent at p = p . 
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Figure 5: Comparison of the "left" and "right" side of the QCD sum rules 
(15) in the vacuum (solid and dashed line) and at normal nuclear density (dot 
dashed and long dashed line) as a function of the Borel mass parameter M. 



Figure 6: The Isobar model for the p-wave interaction of the p- meson (wavy- 
line) with the nucleons (solid line). The double line are the resonances, where 
explicitly include the A(1232), iV*(1720) and A*(1905). 



5 Can one observe chiral restoration? 

We now want to compare our result of the correlation function with the dilep- 
ton data measured by CERES [§. Since the dileptons at this high energy 
collisions carry momenta which are of the same order as their invariant mass, 
we also have to include the p-wave modification of the correlation function. 



Friman and Pirner |14| used an isobar model shown in fig. y in order to take 
the p-wave interaction of the p-meson with the nucleons into account. They 
include in particular the iV*(1720) and the A* (1905) resonance. For the other 
vector mesons no such resonances are known. This gives rise to a momentum 
dependent vector meson-nucleon scattering amplitude Tvn(w, q) which we also 
include in eq. ([3~6|) . Under the assumption that the change of the spectrum is 
mainly due to the density, we can use the resulting ImIT(g, p, T = 0) as in- 
put for eq.([TJ) and determine the thermal dilepton rates. Integrating this over 
the different momenta of the dileptons, the volume V(t) = N%/ p{t) and the 
time evolution of the fireball, we are able to estimate the rate of dilepton pairs 
emitted in a heavy ion collision of sulfur on gold at CERN: 

rf2iW - ^ dtV(t) ■ (20 ) 



dr]dm e + e - Jo J qo d^xd^q 

where A(q) is the acceptance of the CERES detector. Now we just have to 
assume a suitable density and temperature profile which we adapt from the 
usual transport code calculations |TjJ: 

T(t) = (T* - T°°)e- t/Tl + T°° (21) 

with an initial temperature T % = 170 MeV. For the final temperature we take 
T°° = 110 MeV and we use T\ = 8fm. For the density profile we use 

pit) = Pi exp(-t/T 2 ); (22) 

with pi = 2.5p and r 2 = 5fm. We show the result in fig.[7]a and compare it 
with the case where we use the free spectrum at p = 0. One clearly observes 
a strong improvement as compared to the free case. 
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Figure 7: Comparison of the dilepton rates from sulfur on gold collisions 
measured by the CERES collaboration || with 

a) a hadronic system assuming medium or vacuum properties 

b) a system of uncorrelated quark and gluons or a mixed system of this and a 
hadronic phase. 

Now we want to study whether chiral restoration might already occur in 
such collisions and assume, as an extreme limit, as system of free, uncorrelated 
quarks and gluons. Then the spectrum of the correlation function would have 
the constant value of the perturbative plateau in eq . (|i~2"D . We therefore assume 
— 12ir / q 2 lmll(q 2 ) = 5/3 below the strange quark threshold and 2 above. The 
result is shown in fig.[7|b using otherwise the same temperature and baryon 
density profile as in the hadronic case. The consistency with the data is as- 
tonishing. It is clear, however, that in a realistic scenario both hadronic and 
quark-gluon phases must be present. If the initial fireball temperature is as 
large as T % = 170MeV, the system should have already a substantial quark- 
gluon component. We therefore assume next that the uncorrelated quarks 
exist at temperatures above 140 MeV. Below that temperature, the hadronic 
scenario sets in. The result is also shown in fig.|7|b. 

Considering roughly comparable magnitudes of dileptons produced in both, 
the hadronic scenario and the phase of uncorrelated quark and gluons, the real 
problem is how one can distinguish them. It is clear that at the moment 
the precision of the CERES data opens no way to separate them. A possible 
signal might be the suppression of the u- and 0-meson peak in the quark-gluon 
phase. But also a strong broadening of these resonances in the hadronic phase 
can lead to a similar suppression. It is therefore essential first to study the 
hadronic properties more precisely, and under more moderate conditions. 

One way is to look for bound u-meson states. In our hadronic model 



for example we observe a strong attractive potential for the cu-meson in nu- 
clear matter. This behaviour is also supported by QCD sum rules. Such an 
attractive potential would lead to bound ui meson states already in a light 
nucleus like Lithium. Such states could be produced via a transfer reaction 
like d + A — > 3 He + (A — 1) + uj [20|. Here the incoming deuteron d picks up 
one proton from the target nucleus A and is detected as 3 He in forward direc- 
tion. During this process the u;-meson is produced and and possibly bound 
in the residual (A — l)-nucleus. The preferred deuteron beam energy is such 
that the u is produced with minimal possible momentum. The momentum 
and energy conservation determines then the energy of the bound a;-meson. 
A typical spectrum which one expects on the basis of a realistic distorted 
wave Green's functions calculation, is shown in FigM We recognize that the 




Figure 8: Typical spectrum of the transfer reaction d+A — > 3 He+(A — l)+u. 
We consider the results from the free potential (dashed line) and from our in- 
medium potential (solid line) 

increased medium width broadens the single bound uo strongly and makes it 
difficult to identify separate levels. Nevertheless, one would expect to observe 
a strong downward shift of the spectral weight which would be a signal for an 
attractive medium potential. Of course, the differential cross section for this 
process is small, and a detailed examination of competing background (not 
taken into account in Fig. 8) is necessary. 



6 Conclusions 

We have compared the dilepton rates arising from a hadronic model with 
those coming from a phase of uncorrelated quarks. Both scenarios describe 
the CERES-data well and can not be distinguished at the present level of 
accuracy. More precise data under less extreme thermodynamic conditions are 
needed, e.g. from the forthcoming HADES experiments at GSI which permit 



better e + e resolution. One should also study the possibility of bound states 
of cu-mesons in order to gain more information and better understanding. 



References 

[1] T.D. Cohen, R.J. Furnstahl and D.K. Griegel, Phys. Rev. C 
45(1992)1881. 

[2] J. Engels, F. Karsch and K. Redlich, Nucl. Phys. B 435 (1995) 295. 

[3] G. Boyd et al, Phys. Rev Lett. 75(1995) 295. 

[4] M. Lutz, S. Klimt and W. Weise, Nucl. Phys. A 542 (1992) 521. 

[5] G. Agakichiev et al., Phys. Rev. Lett. 75 (1995) 1272. 

[6] M. Masera et al., Nucl. Phys. A 590 (1995) 93c. 

[7] G.E. Brown and M. Rho, Phys. Rev. Lett. 66 (1991) 2720. 

[8] K. Saito and A.W. Thomas, Phys. Rev. C 52 (1995) 2789. 

[9] T. Hatsuda and S.H. Lee, Phys. Rev. C 46 (1992) R34. 

[10] M. Herrmann, B.L. Friman and W. Norenberg, Nucl. Phys. A 560 (1993) 
411. 

[11] G. Chanfray and P. Schuck, Nucl. Phys. A 555 (1993) 329. 

[12] F. Klingl, N. Kaiser and W. Weise, Nucl. Phys. A 624 (1997) 527. 

[13] S. Leupold, W. Peters and U. Mosel, |nucl-th/9708"0T6| . 

[14] B. Friman and H.J. Pirner, Nucl. Phys. A 617 (1997) 496. 

[15] F. Klingl, N. Kaiser and W. Weise, Z. Phys. A 356 (1996) 193. 

[16] CD. Dover, J. Hiifner and R.H. Lemmer, Ann. Phys 66 (1971) 248. 

[17] F. Klingl, T. Waas and W. Weise, ^iep-ph/97092T0| , subm. to Phys. Lett. 
B(1997). 

[18] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 147 
(1979) 385. 

[19] R. Rapp, G. Chanfray and J. Wambach, Nucl. Phys. A 617 (1997) 472. 
[20] A. Gillitzer, R.S. Hayano et al., GSI-Proposal # S214. 



